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Since the computer's inception, the cartographic industry has profited from 
repeated hardware and software advances. Computer-aided mapping systems 
effectively reduce labor costs associated with the production and maintenance of 
cartographic documents. Orbiting satellites widen land coverage to encompass 
the entire globe. Significantly improved remote-sensing systems extend the 
cartographer’s field of sight beyond the visible spectrum. Yet despite such 
progress. the process of cartographic capture and encoding remains a continuing 


concern. 


A. WAT A ENTRY 

Current mapping systems lack the capability to automatically enter and 
classify data from the source. Most employ electromechanical digitizers with 
which human operators manually trace lineations and boundaries. recording (X.Y) 
coordinates at defined intervals or pivot points. Collection of information in this 
format is commonly referred to as vector mode storage. Numeric keys. located on 
hand-held devices. assist feature labeling with specified category codes. The 
above procedure is often a time-consuming, fatiguing and error prone operation. 

Alternatively, manufacturers offer automatic digitizing facilities promoting 


raster mode storage. In raster mode, images are divided into scan lines and 


further subdivided into small picture cells. Information is stored in a sequential 
fashion by horizontal scan line pattern arranged from the top of a picture 
downwards. If retention of additional area characteristics is desired. manual 
intervention, tiring by nature, is necessary to separate and identify individual 
components. 

The market also supports automatic line-following digitizers capable of 
sensing and tracking linear features with servomechanisms. Aimed at alleviating 
the attentive and tedious handling common to other methods. these systems have 
limited use and require high quality maps as input. Uncertainties arising at 
intersecting features are resolved through the assistance of an operator. : 

The solution to automatic data input and feature identification from paper 
maps presents a most difficult challenge. The production of a fully automated 
system may possibly never be attained. Present day manufacturers have already 
acceptably confined their systems’ capabilities to only certain sub-tasks. Two 


major reasons account for the limited progress seen in the area of data capture: 


sx the market is characteristically small and undercapitalized. offering little 
growth potential. 


mz hardware facilities have lacked sufficient computing power and display 
resolution capabilities to handle the demanding load of processing and 
manipulating cartographic data. 


B. MOTIWATION 

Many government mapping agencies and commercial map companies have 
realized the potential benefits of maintaining digital cartographic data bases. The 
U. S. Geological Survey. for example, is in the process of converting all of its 
United States quadrangle series maps into digital form. Topographic data is being 
accumulated into separate categories such as digital elevation, surface 
hydrography, public land-survey network, geographic names and other classes. 
Specialized projects, requiring digitized map files. are newly emerging. For 
instance. work conducted by Richbourg, Rowe and Zyda [Ref. 1] on solution 
techniques for two-dimensional route planning problems for mobile autonomous 
vehicles identifies a need for topographic information grouped by terrain speed 
regions. 

The inability of current data bases and mapping systems to fulfill the latter 
examples requirements served as a catalytic force in reopening the study for 
improving computer-assisted cartographic entry and encoding techniques. A 
discussion of the processing system under development is presented in Chapter 
Two. The major goal of the research is to offer an alternative approach from the 
popular interactive digitizing methods used in handling data entry from published 
maps. Although it is difficult to design an accurate, reliable and low cost scheme. 


the proposed approach seems promising. 


C. SCOPE OF fii si bar 

This study develops and examines an interactive editing system for 
polygonally encoded data bases generated from the cartographic processing system 
described in |Ref. 2]. The editor, serving as the final step in the overall project, 
provides the user with the capability to update and modify topographic 
information. A variety of common processing and digitizer induced errors. 
including misencoded coordinates, line noise and superfluous points can also be 
corrected. The editor's algorithms were implemented using the C programming 
language with calls made to available Graphics Library routines. A Silicon 
Graphics. Inc. IRIS (Integrated Raster Imaging System) Turbo 2400 workstation 
was chosen for this study as it supports the high resolution and fast transmission 


rates needed for manipulating cartographic data. 


D. ORGANIZATION 

Chapter Two reviews the complete cartographic data processing pipeline 
under development at the Naval Postgraduate School. Each step within the 
overall package is briefly discussed. Chapter Three reports on the available 
hardware and software research facilities utilized. Chapter Four introduces a 
graphical editing system. the topic of this study. Chapter Five addresses the 
techniques employed for managing user revisions. The algorithms used for 
transforming an updated map image from the screen into stored sie are also 


covered. Finally Chapter Six is an examination of the presented interface’s 
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limitations and weaknesses with suggestions for future enhancements and 
improvements. Attached as appendices are before and after snap shots taken from 
sample editing sessions. a list of USGS map symbols supported by the editor and 


the algorithms utilized for updating a cartographic data base. 
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TI. A CARTOGRAPHIC DATA PROCESSING SYSTEM 


A. OVERVIEW 

The integrated cartographic processing package, developed at the Naval 
Postgraduate School. is organized as a set of five stages, each performing a 
distinct function. <A collection of methods and techniques. well-known in some 
cases and novel in others. assist in making this prototypical model operational. 
The component steps are shown schematically in Figure 2.1 in the order they are 
generally executed. Two of the stages (items 1 and 2) deal primarily with the 
selection and preparation of maps for use in subsequent steps; two (items 3 and 4) 
handle actual digital information processing and the creation of a computerized 
data base: and the remaining one (item 5) conducts ein and editing of files 
stored on disk. Steps three and four are the subject of research conducted by 
Diehl {Ref. 2] while step five is the focus of this study. Each stage of the project 


is briefly described below. 


B. PRE-PRINTED MAP SELECTION 

Prior to implementing any cartographic information processing system. 
selection of a standardized mapping scheme is necessary. For this research. 
primary input consisted of sectional pieces from U. S. Geological Survey (USGS) 


Quadrangle (7.5-minute x 7.5-minute) series maps with conventional units at a 
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1:24.000 scale. Widely available and accepted. each USGS map complies with 
established specifications in regards to size. scale. symbolism and content |Refs. 3 


and 4]. The use of color assists in discriminating the various cartographic features 


(see Table 2.1 below). 


TABLE 2.1 
USGS MAP COLOR USAGE SUMMARY 


Application Area 


BoA Cultural features such as roads and 
buildings. 
BLUE Hydrographic features such as lakes. 
a. rivers, wetlands and reservoirs. 


BROW! Hypsographic features shown by contour 
ee ee 
GREEN Surface cover including woodlands, scrub. 
POMS | asartne eats 
PURPLE Features added from aerial photographs 
during map revision. The changes are 
not field checked. 


Important roads and public land survey 
system. 


WHITE Nonvegetative features such as barren 
waste areas and piles. 





C. DIGITIZATION 

The raw pictorial representation of the paper map is converted into 
computer-readable form for subsequent processing through the use of a video 
digitizer camera. Digitization comprises the scanning of a document to resolve 


line and character information into small picture elements, or pixels. The optical 
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response at each pixel controls detection circuitry for the generation of 


positionally defined digital output. 


An EYECOM! Picture Digitizer and Display unit interfaced to a Digital 
Equipment Corporation (DEC) PDP-11 computer was utilized for this project. 
Additional auxiliary lighting and filters (red, green and blue) complement the 
setup. Files produced after digitization are downloaded onto a DEC VAX 11/780 
computer. The File Transfer Protocol (FTP) program on the connecting 
ETHERNET allows file movement to the final destination host. a Silicon 


Graphics, Inc. IRIS Turbo 2400 workstation. 


D. IMAGE PROCESSING 
1. General 

After ees inine a Maps contents into digital form, ideally one applies 
an image processing system which promotes automatic component recognition. A 
component is any feature that is judged to be significant at a chosen scale and 
that exists within an area. As discussed above, present day systems fall short in 
reaching this optimum capability. Most if not all commercially available systems 
require varying degrees of human intervention to demarcate cartographic 
information. Although the market supports some automatic line-following 
systems, these tools are useful only for certain types of features and only on maps 


with high qualitv line work. In general, human assistance is needed to bridge 


1 A trademark of Spatial Data Systems, Inc. 
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interruptions. to resolve uncertainties at intersections marked by unreliable 
directional trends and to affix identifying codes for individual features. 

A technique proposed by Diehl! [Ref. 2| attempts to resolve this noted 
deficiency. Although in a prototypical state, his method offers an alternative 
avenue for further research exploration in solving this most dificult problem. An 
overview of his research follows. 

2. The Diehl Feature Extraction System 

From a black box viewpoint (see Figure 2.2). Diehl’s system analyzes 
individual pixels of an input video image and through heuristic reasoning, 
categorizes each pixel into a subset of colors synonymous to those comprising the 


initial paper map. In the digitization process. each pixel has the potential of 


adopting one of over sixteen million color identities” under the RGB color system 
- 1. e. 24 bits per pixel. with 8 bits each for red, blue and green intensity values. 
Determination of actual pixel color assignment is derived by sample readings 
taken eel contiguous map areas. This averaging procedure accounts for 
color variations observed in pixels of seemingly uniform regions. Points lying in 
oe proximity to boundaries separating dissimilar colors often assume fuzzy values 
while interior ones approximate original map colors to within a measurable 
deviation. Preliminary tests show Diehl’s work to be highly accurate in correctly 


assigning pixels of internal feature sections with small error 


* All colors in the spectrum can be formed by blending the additive primitives (red, green and 
blue). 
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rates occurring near obscure border lines. The resulting data base generated from 
this processing scheme can either be retained as is or be funneled through the 


Data Compression Stage to optimize storage space. 


E. DATA COMPRESSION 
Once the map has been segmented into one of seven equivalence classes, an 
appropriate storage format (raster versus vector) must be chosen. Several factors 


were considered during the selection process : 
sez speed tolerance for terminal redisplay and screen refresh. 
sx degree of resolution needed to maintain accuracy, 
re retrieval rate for information of interest, 
ce ease in data manipulation during editing sessions, 
me storage resource requirements, 
ce quantity of information to be retained per unit of measure, and 


re applicability of format to problem area. 


Since our interest centers on the generation of data bases for use in autonomous 
vehicle route planning, vector storage proved to be most economical. 

A triangular tessellation. formulated by Diehl [Ref. 2]. serves as the 
fundamental data storage mechanism. In general, the process consists of two 


repeatedly called parts : 
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c= identification of the next homogeneous area. and 


cz division of the identified region into polygons. 


A detailed discussion of the algorithms involved is presented in Chapter Five. 
These algorithms were slightly modified and incorporated into the editing 


system's facility for updating a cartographic data base. 


F. THE EDITOR 

Finally, an editing system is introduced to ensure proper maintenance of the 
data base generated from the previous utility steps. Typically most printed maps 
on the market are several years si and are rarely updated. A computer-assisted 
map data base editor allows easy modification and correction of dated 
topographic characteristics. The digital map data can be altered either to 
eliminate or remedy errors and blunders or to revise established files when new 
updates arrive. Often proved to be time consuming for the user, this is an 
essential step in the overall process of assuring a maps usefulness and accuracy. 
Possible adverse effects could result based on outdated or erroneous information 
portrayed by an unkept data base. The contents of the remaining chapters focus 


on the editing system s design. 
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Il... RESHAR @ ier Ae ile 


A. IRIS WORKSTATION 

The IRIS Turbo 2400 Graphics Workstation, manufactured by Silicon 
Graphics, Inc.. served as the primary tool for this research. By incorporating 
aa built VLSI circuits into its design. the IRIS offers an attractive alternative 
to the more conventional workstations [Ref. 5]. The special-purpose hardware, 
designed to replace less efficient software. yields high processing speeds and 
increased performance reliability necessary for manipulating cartographic data. 
The system notably combines real-time color graphics with Unix operating system 
software and Ethernet network communications. A very high resolution color 
monitor provides crisp displays and fine delineations for even the most complex 


work. Figure 3.1 outlines the distinctive features of the system. 


B. IRIS GRAPHICS LIBRARY 

The IRIS Graphics Library is an extensive collection of utility and graphics 
subroutines providing high level access to the hardware enabling graphical objects 
to be easily manipulated as geometrical objects (points. lines, polygons, etc) 
rather than pixels. A series of coordinate systems and mapping instructions 


provides the user with the capability to define such objects in world space. 
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A subset of the Library's routines were used for this research. They can be 


generalized into the following nine categories : 
- Global state commands initialize the hardware and control global! state variables. 


- Primitive drawing commands draw points, lines, polygons, circles, arcs, and text strings on 
the screen. 


- Drawing attribute commands select characteristics for drawing lines, filling polygons, and 
writing text strings. 


- Coordinate transformation commands perform manipulations on coordinate systems, 
including mapping user-defined coordinate systems to screen coordinate systems. 


- Display mode and color commands determine how the bitplane image memory is used and 
how objects are colored on the screen. 


. Input /output commands initialize and read input/output devices. 


- Object creation and editing commands provide the means to create hierarchical structure 
of graphics commands. 


- Picking and selecting commands identify the commands that draw to a specified area of 
the screen. 


- Teztport commands allocate an area of the screen for writing text.° 


A detailed account on command usage for each of the above classes is covered in 


the IRIS User’s Guide. 





“Silicon Graphics Inc.. JRIS User’s Guide Version 2.1. Page 1-2. 
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Ethernet connection between IRIS’s and VAX’s 





32-bit Motorola 68020 Processor 

1024 X 768 X 32-bit Display Memory 

2 MB of CPU Memory 

144 MB of Disk Storage 

Floating Point Accelerator 

Geometry Pipeline with Geometry Engines and Geometry Accelerators 
6O Hz Non-interlaced Display 

16-bit Z-buffer for Hidden Surface Elimination 
Hardware Smooth Shading 

Unix System V 

IRIS Graphics Library 

Cartridge Tape Unit 

Ethernet to VAX’s 


Digitizer Tablet 


Features of the IRIS Turbo 2400 Graphics Workstation 


Figure ofa 
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he lie EDITING SYSTEM 


AemeGit LING St ARTisD 
The editing system is invoked by a single command to the IRIS workstation : 
% MAPEDIT 
Upon execution. an initial title page is displaved for a few seconds prior to the 
appearance of the main program. Prompts and instructions guide the user 


through the entire operation of the editor. 


B. SCREEN LAYOUT 

The screen is partitioned into five non-overlapping windows (see Figure 4.1). 
Each is reserved for a particular use. Window One serves as the primary work 
bench for area viewing and editing of a sectional map. Window Two exhibits a 
dynamic picture display of the current map image in miniature form. Window 
Three contains any menu options available during a particular point in time. The 
use of pop up menus. described below, assists in reducing screen clutter and user 
distraction. Window Four is restricted to point viewing and editing. A small 
region (20 pixels by 20 pixels) surrounding the desired point is magnified by a 
factor of ten. Finally, Window Five displays pertinent system status information 


and helpful instructions to the user. 
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The editor also emplovs temporary and overlapping windows to allow more 
surface area from a fixed size screen. Designed to preserve partially obscured 
terminal contents. the windows enable viewing of information on drawing 
selections. feature classifications and directory listings. These overlays are 


displaved only upon request. 


© INPUT DEVICES 

User responses are entered either through the IRIS keyboard or three button 
mouse. The keyboard’s function is limited to the entry of input/output filenames. 
desired directories for display. and feature categorization codes. Remaining input 
is regulated by the mouse. 

The three mouse buttons (Left, Middle and Right) offer flexible and easy 
control. By properly positioning the cursor on the monitor screen and 
simultaneously pressing the appropriate button, the user can effortlessly perform 
such actions as the selection of menu choices. areas to be updated and different 
drawing options. Movement of the mouse across a hard surtace dictates the 


cursor s actual location on the screen. 


D. MBNUTSTRUCTURE 
The editing system consists of a series of pop up menus organized in a 
hierarchical structure (see Figure 4.2). Source code for these menus was derived 


by modifying a general purpose menu package provided by Silicon Graphics. 
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The pop up menu routine is very simple. When invoked. a menu is 
automatically drawn on the screen displaying the possible program selections. By 
moving the cursor with the mouse device across an option. that entry is 


highlighted. Upon the depression of any mouse button. three actions occur : 
pa the entry under the cursor at that time is selected. 
re the menu disappears from the screen. and 


r= the program branches to the procedure specified by the option. 


If the cursor lies completely outside the menu field when the button is pressed. 
the menu remains on the screen and the routine continues to wait until an 
appropriate response is entered. 

A description of the operations and capabilities for each of the available 
commands follows. 

1. Area Edit 

This procedure permits large region editing of the map by allowing the 
user to add. delete. move and change displayed components. Refer to the section 
"Point Edit" for modifying individual pixels. 

Features are designed by combining different drawing commands and 
attributes. Other modifications such as deleting, reproducing, moving and 
changing attribute settings require the desired figure be identified. Selection of an 
individual object is possible by positioning the cursor over it and then depressing 


the middle mouse button. The editor will flash the chosen figure and request 


27 


confirmation of its correctness. A rejection by the user causes the system to 
search for the next best solution. This process of searching and requesting 
approval continues until the user is satisfied with an object or when no further 
figures are Pawan If a selection is made. the current editing action is then 
performed. A warning is displayed when no figure is found. 


Six choices are available in the Area Edit Menu : 


(1) Drawing Commands - allow feature sketching using one or more of 


the following shapes’ : 


Circle 
Rectangle 
-- Polygon 
Horizontal Line 
Vertical Line 
Single Line 
Continuous Line 


The last produced shape can either be erased or reproduced without 
entering the "Delete Figure" or "Reproduce Figure" choices of the Area 
Edit Menu. The new shape constructed will be situated at the cursors 
current location for further repositioning on the map. 


(2) Attributes - enable viewing or redesignation of the current setting or a 
specific figures characteristics. This option comprises four subcommands: 


View Attribute Set 
Change Attribute Set 
View Figure Attributes 
Change Figure Attributes 


Selection of the first subcommand provides a display of the current 
attribute settings without modification. 


* CAUTION : The IRIS system does not support concave polygons; drawing these figures 
may produce unpredicatable results. 
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The second subcommand allows any of the current attributes to be 
changed. When the editor is first entered. they are initially set to the 
following values : 


Color - White 

Linestyle - Continuous 
Linewidth - QOne pixel wide 

Fill Pattern - Solid 

Feature Code - 001 (nonvegetative) 


The first four settings are modified by choosing a desired alternative from 
among the visually displayed pop up pallets. There are twelve colors 
(black, <light and dark> blue, <light and dark> brown, <light and 
dark> green. <light and dark> purple. red, yellow (unknown) and white). 
fourteen linestyles. five linewidths and ten fill patterns. Feature codes are 
altered by entering in the new value after a system prompt. Depression of 
the PF1 key displays the entire listing of descriptions with corresponding 


_ codes for review. 


The third subcommand allows one to view a particular figure’s attributes 
without modification. 


To change an individual figure’s attributes. the fourth subcommand is 
employed. After the figure is located. current settings for the figure are 


displayed. Changes are entered in a similar manner as described in the 
process for altering current attributes. 


Delete Figure - removes the selected object from the display. The 
vacancy is filled with the color white and identified as nonvegetative. 


Reproduce Figure - produces a carbon copy of the selected object with 
identical attribute settings. 


Move Figure - repositions the selected object to any desired region on 
the map. The cancel option returns the object to its original location. 


Change View - refer to the section "Change View" discussed below. 
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Point Edit 


Area Editing commands can become quite cumbersome when performing 


slight alterations to existing features or when identifying unknown regions around 


border lines. The Point Edit facility offers a viable alternative. With this 


procedure, one can easily modify individual pixels within a map display using the 


following two commands : 


(1) 
(2) 


Point Zoom-In - refer to the section "Change View". 


Attributes - allow viewing or redesignation of the color and/or feature 
description of the current setting or of any specified pixel within the map. 
The desired point to be updated is first located with the Point Zoom-In 


. command. Selected pixels assume the currently specified color and feature 


categories. By default, these settings are initialized to nonvegetative 
(color = white. feature code = 001). The editing color may be changed at 
any time from a pallet of map colors displayable upon request. Feature 
description codes are modified by entering in the new value after the 
system prompt. The entire listing of description codes can be viewed by 
depressing the PF1 key. Changes recorded on the magnified area are also 
reflected in the current map view. 


Change View 


Magnification of different map sections may be necessary to assist the 


human eye in correctly perceiving and locating small features requiring updating. 


Three options are provided by the editor : 


(1) 


Area Zoom-In - 2x and 4x magnification of selected user-defined areas. 
The current map view is overwritten by the enlarged area. (screen 
location : Window One). 


Point Zoom-In - 10x magnification on a 20 pixel x 20 pixel area 
surrounding the desired point. (screen location : Window Four). 
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(3) Reset - restores the screen to the original map view. 


Exiting this procedure. when called from the Main Menu. invokes the Reset 
command. 
4. Read ina File 
This choice allows entry of previously stored map data for further 
editing. After display of the system prompt. the user types in the desired 
filename at the keyboard. If the file exists. the data is drawn on the screen 
overwritting any prior map images; otherwise, an appropriate error message 
appears with control returning to the Main Menu. The user may cancel any 
inadvertent calls made to this routine by using the Escape key on the keyboard. 
A directory listing facility is also available. By depressing the PF1 key. 
a prompt is shown requesting the desired directory name; the user's current 
directory serves as default. This facility is capable of listing all files in any 
directory not read protected. 
5. Write toa File 
To save changes entered by the user. a method is necessary to capture 
the screen display as stored data within a file. This command selection provides 
such a capability. The projected map image is transformed into similarly colored 
areas which are then subdivided into polygons before being written to the created 


file. Specification of an existing filename results in the destruction and 
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overwritting of any old contents. The actual algorithm used for this procedure is 
presented in detail in Chapter Five. 

The process of designating an output filename in addition to the Escape 
and PF1 key functions are identical to those previously described (refer to the 
section "Read in a File"). 

6. Exit System 

Selection of this option terminates the editing session. Changes not 
previously saved are lost forever. The graphics hardware system is returned to its 
initial state. Input devices are unqueued. the viewport and textport are restored 


to full-screen. and the communication buffers are flushed. 


BE. SYSteneen OME: 

System information and prompts appear on the screen at various times 
during the execution of the editor. The prompts serve as a guide through an 
editing session without reference to an operating manual. These helpful messages 
reduce the amount of information the user must memorize while keeping him 


informed of what is happening now and what to expect next. They include : 


(1} Menus - allow selection among multiple program alternatives. (screen 
location : Window Three). 


(2} Key Function Instructions - display the active keyboard keys at any 
given time. The keys are identified by dark blue squares with an 
abbreviation of the key inside (for example. ESC signifies the Escape 
Key}. Keys not listed are inactive; pressing them has no effect. (screen 
location : Window Five). | 
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(3) Mouse Function Instructions - list the outcome of depressing anv of 
the three mouse buttons. The buttons are represented by dark blue 
squares with a letter ( L for left M for middle and R for right ) inside. 
Buttons not shown are inactive: pressing them has no effect. (screen 
location : Window Five). 


(4) System Status - report on the program’s reaction under the following 
conditions : reading in a file. writing out a file. delay due to processing. 
and no action taken. (screen location : Window Five). 


(5) Error Messages - warn when erroneous or adverse conditions are 
present (for example. invalid input). Depending on the situation, the 
system tries to recover by returning to its last working state. (screen 
location : Window One). 


(6) Other Prompts - request additional input information from the user 
such as a filename. directory name and/or feature category code. (screen 
_ location : Window One). 


fete WR SUPPORTED 

Cartographic components are designed and updated using ine Area and 
Point Edit options described above. The drawing commands and attributes, 
provided by the editing system. enable the user to generate most of the symbology 
related to USGS maps (see Appendix B). 

A three digit description code assists in feature identification. During the 
pre-processing stages. all components are initially classified into one of the 


following broad categories based on color: 


Color Code Description 
Yellow 000 Unknown 

White 001 Nonvegetative 
Black 002 Man made object 
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Red 003 Road 


Blue 004 Body of Water 

Brown 005 Contour line 

Green O06 Vegetation 

Purple 007 Feature from aerial photo 


These initial settings can be redefined with either edit facility. The remaining 
possible feature descriptions with corresponding program defined codes are listed 


in Appendix B. 


G. FILE FORMAT 

The cartographic data base used for this study was generated from the 
processing scheme presented in Chapter Two. The data within each file represents 
a section from a USGS map stored in binary format. 

Each map section is divided into a number of uniformly colored areas. An 
area is identified bv a set of short integer codes representing its color. description 
and allowable speed zone. Each area is further subdivided into ; number of 
polygons as determined by the Data Compression Algorithm. The polygons are 
defined by short integer codes representing their type (point. line. polygon, filled 
polygon) and vertex locations. An end of area marker separates each successively 
defined area region. Figure 4.3 depicts a typical file format. 

The editor extends the seven map colors to twelve. Light and dark colors for 
green. blue, purple and brown support cartographic feature contrast in given 
regions and have no significance with respect to classification. Yellow identifies 


unknown areas. 
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@ 
color description speed 


draw command vertices 


draw command Vertices 


end of area marker 





=. 
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* 
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@ 
S 
@ 
end of file 


color ai: 
HOt ot draw commands verticies 


black 


‘ point xlyl 
white 


line xlyl x2y2 


red 
xlyl x2y2 x3y3 


polygon 


light green 
filled polygon xlyl x2y2 x3y3 


dark green 
light blue 
dark blue end of area 


light purple ~1 (not used) 


speed 


dark purple 


light brown description codes 


(see appendix B) 


dark brown 


O 
1 
2 
3 
4: 
5 
6 
7 
8 
9 
10 
11 


yellow (unknown) 





Figure 4.3 Typical File Format with Key 
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Ve MAWWAGING USER REVISIONS 


A major concern in the design of any man-machine interface is the ability to 
insure real-time response rates. Users by nature become very anxious when 
reactions to their requests do not occur within a few seconds of input. Even 
prompting messages, indicating processing is in progress. lose their effectiveness 
after a period of time. The user may begin to fidget and start pressing buttons 
and "trying stuff", perhaps leading to disastrous results. In order to minimize 
such processing delays experienced during an editing session, an internal indexing 


scheme was incorporated into the editors design. This plan. in addition to the 


method used for updating a data base. are presented below. 


A. INTERNAL INDEXING SCHEME 
As mentioned in Chapter Four. each map is divided into several uniformly 


colored areas which are further subdivided into a number of polygons. Each area 


is defined as an object.° Due to the feature complexity of most cartographic 
documents. it is not unusual to have in excess of 5000 objects defined at any point 
in time. This large number places a strain on the IRIS'ss memory management for 


locating a particular object quickly. Response time slowly decreases at an inverse 


> An object, the graphics packaging mechanism provided on the IRIS, is composed of a 
sequence of graphics commands. 


rate to the number of objects created. To help alleviate this problem. an artificial 
indexing scheme (see Figure 5.1) was added to the editing system. 

The scheme consists of an indexing array subscripted from O through 3277. 
As an area of the map is entered into memory. it is assigned a short integer 
identification number. The area is then placed into a linked list, arranged in 
descending order. at an array location equal to the absolute value of the index 
number divided by ten. 

To locate an object for updating or modifying, the pick mechanism from the 
IRIS graphics library is employed. The picking facility allows the user to identify 
an object by maneuvering the cursor as a pointing device. Objects intersecting a 
small rectangular area around the cursor are considered "hits" and their 
identification numbers. assigned during input. are placed into a names stack. The 
indexing scheme. which reduces search time to a maximum of 20 objects. can then 


be used to quickly locate each object by number for manipulation. 


Pott OS tiNG THE DATABASE 

Upon completion of an editing session, changes to a map are saved by using 
the "Write to a File" option from the Main Menu. The four algorithms, listed in 
Appendix C. comprising this command selection include : Update Database. 
Data Compression, Expand Area and Subset Area. The latter three were 


written and used by Diehl in the Data Compression phase of the cartographic 
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processing package. Slight modifications were added to enable handling of a 
larger set of colors. Each of the algorithms are described below. 
1. Update Database 

Since the primary use for these data bases is in autonomous vehicle route 
planning, it is desirable for components to be uniquely represented in polygonal 
form. Duplicated areas. as seen in the case of overlying features. need to be 
eliminated. The Update Database Algorithm provides this capability. 

The program itself assigns a set of unique color shades. representing the 
twelve possible map colors. to each feature description. During an editing 
session, default colors are actually used for display and to reflect entered changes. 
This algorithm replaces the default color for each area of the map with an 
appropriate color shade as dictated by the area’s description code (see Figure 
5.2a). After reassignment. the map is redrawn. The new map image is 
subsequently read off the screen and funneled through the Data Compression 
Algorithm (see Figure 5.2b). 

2. Data Compression 

The entire map surface is constructed as a grid, each cell being one pixel 
wide. A duplicate setup. hereafter known as the scratch pad. has all its cells 
initially set to "O". Working in a sweeping motion from left to right and top to 
bottom on the proposed map layout, the algorithm checks for unexamined pixels. 
When one is found, the process converts the pixel’s color shade back to the 


default color (see Figure 5.2b) and writes the color. description code and speed 
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code to a file. The Expand Area and Subset Area algorithms are invoked prior to 
outputting an end of area marker. This cycle persists until each pixel of the map 
has been surveyed. 
3. Expand Area 
This algorithm focuses on a three by three subgrid successively 
encompassing each cell of the map. Three actions occur when any of the 
surrounding eight neighbors have the same color as the center. er each similarly 


colored pixel : 
sa the pixel is tagged as used to avoid duplicate inspections, 


mr a flag is set to "1" in the appropriately corresponding position of the 
scratch pad, and 


se the pixel’s location is placed in a queue for later examination of its 
neighbors. 
The interrogation process continues until the queue is empty. resulting in a 


uniformly colored region marked by set flags in the scratch pad. 


4. Subset Area 
The identified homogeneous area marked in the scratch pad is then 
subdivided using the ten possible geometric shapes shown in Figure 5.3. Items a 
through e represent right equilateral triangles. items f through i describe lines, and 
item j defines a point. |Ref. 2] addresses the actual selection criteria for these 


forms and the process involved in dividing intact area regions. Upon completion 
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of the subdivision. scratch pad flags are reset to "O" and the area identification 
phase is reexecuted. 

The above method is a simple scheme allowing compact storage for 
processed and edited information. The resulting data base itself contains a 
cartographic image divided into homogeneous areas which are further subdivided 
into polygons. Each polygonal side and/or vertex is uniquely defined avoiding 
redundant representations. One and only one feature classification is attached to 


each polygon. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

This study examines an interactive editing system for maintaining data bases 
generated from a cartographic processing system currently under development at 
the Naval Postgraduate School. The amount of updating necessary on any stored 
map depends on the care taken in digitization. the accuracy of the image 
processing phase and the level of quality required in the finished product. In 
general. the time devoted towards this effort is proportional to the maps 
complexity. If each step involved in the generation of the data base has been 
carefully executed and there are no mistakes in content. the editing is mostly 
cosmetic. 

The editor was written in the C programming language and contains 
approximately 10.000 lines of code. Driven by a set of pop up menus arranged in 
a hierarchical structure. the system has been designed to handle modifications to 
either large map areas or individual points. The latter capability is useful! in 
correcting processing and digitizer induced errors around linear features. A mouse 
served as the primary input device and proved to be fairly accurate in 
manipulating and designing cartographic data. Display colors approximate those 


used on the paper map. 


Preliminary tests indicate the editor's performance depends on the maps 
complexity and the hardware’s internal memory size. The system performs best 


using an IRIS with six megabytes of memory and a floating-point accelerator. 


B. LIMITATIONS AND RECOMMENDED EXTENSIONS 

The editor developed in this study merely scratches the surface of the 
problem in applying interactive techniques to update and modify cartographic 
data bases. Several limitations prevent the present system from operating at an 
as yet practical level in a real time environment. A few of these limiting factors 


with possible research extensions are discussed below : 


1. Each preprocessed map contained within the database defines a 280 x 440 
pixel area. This self-imposed size constraint is addressed in [Ref. 2]. Many 
cartographic features commonly span across several adjoining map 
sections. Ideally, any portion of these features should be readily available 
for editing. reflecting all updates. The current editing system is unable to 
handle changes occurring over several neighboring regions. Modifications 
are instead restricted to within one defined map section at any given time. 


EXTENSION : Addition of a facility to merge adjacent sectional map 
pleces and modification of the editor to handle the larger areas. 


to 


USGS cartographic features served as a standard for this research. They 
are designed by combining different colors. linewidths, linestyles and fill 
patterns. The editing system. however. lacks appropriate checks for 
preventing entry of illegal combinations which may render a data base as 
unusable. 


EXTENSION : Incorporation of proper error checking to the editor to 
minimize such mishaps. 


Some USGS cartographic features are not supported by the editor due to 
their inherit complexity and the difficulty of approximating their display 
by setting bit maps such as with Tailings. Intricate Surface Areas and 


co 


46 


Gravel Beaches. Others require several overlays in order to capture the 
proper design. The latter case. in particular. is very tiring to the user and 
ineficient for the computer. Additionally. the algorithm employed in 
updating the database (discussed in Chapter Five) fails to preserve the 
integrity of areas composed of complex fill patterns. Instead. these areas 
are needlessly splintered into smaller regions. Many of the standard 
cartographic symbols look nice on paper but significantly tax memory 
requirements in duplication efforts on the computer. 


EXTENSION : Development of a more efficient graphical representation 
for mapping symbols. The proposed scheme should allow one to depict 
the various natural. man-made and aerial features. The result achieved 
should be equal to traditional cartographic techniques as far as accuracy. 
integrity of design and quality of reproduction. 


4. The Data Compression Algorithm used for subdividing a map produces 
areas composed of many irregularly sized polygons having no common 
sides. Due to these variations. editing of sub-areas tends to be awkward. 


EXTENSION : Modification of the Data Compression Algorithm to 
produce more equally sized subdivisions. 


5. Curved features (contours, winding roads/rivers, etc) can only be 
approximated by using several straight line segments. The accuracy of 
these features is dependent on the users drawing skills. 

EXTENSION : Development of a procedure capable of generating 


accurate curved lines without excessive memory or execution 
requirements. 


6. No text labeling procedures are provided by the editor. 


EXTENSION : This facility is not necessary unless the editor’s role is 
extended to uses requiring that capability. 


Despite a number of limitations. the editing system illustrates the useful role 
interactive computer graphics plays in maintaining geographic data bases. 


Additional studies are needed to expand and improve upon the current system s 
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capabilities of allowing new editions and _ revisions without’ expensive 


recompilation. 
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APPENDIX A — SAMPLE EDITING SESSION 
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APPENDIX B — USGS TOPOGRAPHIC MAP SYMBOLS 


[Ref. 6) 

CODE DESCEIE TIO 

010 CONTROL DATA AND MONUMENTS 

O11 Aerial photograph roll and frame number t 

012 Horizontal control: 

013 Third order or better, permanent mark z, 

014 With third order or better elevation a 

O15 Checked spot elevation | 

016 Coincident with section corner tr 

O17 Unmonumented i 

018 Vertical control: 

019 Third order or better. with tablet 

020 Third order or better. recoverable mark + 

021 Bench mark at found section corner 7 

22 Spot elevation + 
- 023 Boundary monument: 

024 With tablet a5 

025 Without tablet F 

026 With number and elevation , 

027 U.S. mineral or location monument 

030 BOUNDARIES 

031 National 

032 State or territorial 

033 County or equivalent 

034 Civil township or equivalent 

035 Incorporated city or equivalent 

036 Park. reservation. or monument 

037 Small park 

040 LAND SURVEY SYSTEMS 

041 U.S. Public Land Survey System : 

042 Township or range line 

043 Location doubtful 


O93 


CODE 


044 
045 
046 
047 
048 
049 
050 
O51 
052 
053 
054 
055 


O80 
O81 
082 
083 
O84 
O85 
O86 


DESC Tie; 


Section line 
Location doubtful 
Found section corner 
Found closing corner 
Witness corner 
Meander corner 
Other land surveys: 
Township or range line 
Section line 
Land grant or mining claim 
Monument 
Fence line 


ROADS AND RELATED FEATURES 
Primary highway 
Secondary highway 
Light duty road 
Unimproved road 
Trail 
Dual highway 
Dual highway with median strip 
Road under construction 
Underpass 
Overpass 
Bridge 
Drawbridge 
Tunnel 


BUILDINGS AND RELATED FEATURES 
Dwelling or place of employment : small 
Dwelling or place of employment : large 
School! 
Church 
Barn. warehouse. etc.: small 
Barn, warehouse. etc.: large 


o4 


CODE 


O87 
O&8 
O89 
O90 
O91 
092 
093 
094 
O95 
096 
O97 
O98 
O99 
100 
POL 
ioe 
103 


110 
el 
ee 
ele 
114 
11s 
BIG 
AT 
is 
iat 
120 


+ 
oo 
@) 


_— 
Co Co Co Co CO 
Co 


_ _ 
Ho 


— 
on 


DESCRIPTION 


House omission tint 
Racetrack 

Airport 

Landing strip 

Well (other than water) 
Windmill 

Water tank: small 
Water tank: large 
Other tank: small 
Other tank: large 
Covered reservoir 
Gaging station 
Landmark object 
Campground 
Picnic area 
Cemetery: small 
Cemetery: large 


RAILROADS AND RELATED FEATURES 
Standard gauge single track 
Station 
Standard gauge multiple track 
Abandoned 
Under construction 
Narrow gauge single track 
Narrow gauge multiple track 
Railroad in street 
Juxtaposition 
Roundhouse and turntable 


TRANSMISSION LINES AND PIPELINES 


Power transmission line: pole 
Power transmission line: tower 
Telephone or telegraph line 
Aboveground oil or gas pipeline 
Underground oil or gas pipeline 
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CODE DESC Re, 


140 CONTOURS 

141 Topographic: 

142 Intermediate 

143 Index 

144 Supplementary 

145 Depression 

146 Cut 

147 Fill 

148 Bathymetric: 

149 Intermediate 

150 Index 

151 Primary 

ho Index Primary 

153 Supplementary 

160 MINES AND CAVES 
161 Quarry or open pit mine 

162 Gravel. sand. clay. or borrow pit 

163 Mine tunnel or cave entrance 

164 Prospect 

165 Mine shaft 

166 Mine dump 

167 Tailings 

170 SURFACE FEATURES 
ge. Levee 

172 Sand or mud area. dunes. or shifting sand 
173 Intricate surface area 

174 Gravel beach or glacial moraine 

io Tailings pond 
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CODE 


180 
181 
182 
183 
184 
185 


190 
191 
192 
193 
194 
195 
.. 196 


200 
201 
202 
Zo 
204 
205 
206 
Zs 
208 
209 


210 
Zit 
Zl 
213 
214 
215 
216 


DESCRIPTION 
VEGETATION 
Woods 
Scrub 
Orehard 
Vinevard 
Mangrove 


MARINE SHORELINE 
Topographic maps: 
Approximate mean high water 
Indefinite or unsurveyed 
Topographic-bathymetric maps: 
Mean high water 
Apparent (edge of vegetation) 


COASTAL FEATURES 
Foreshore flat 
Rock or coral reef 
Rock bare or awash 
Group of rocks bare or awash 
Exposed wreck 
Depth curve 
Sounding 
Breakwater, pier, jetty. or wharf 
Seawall 


BATHYMETRIC FEATURES 
Area exposed at mean low tide 
Sounding datum 
Channel 
Offshore oil or gas: well 
Offshore oil or gas: platform 
Sunken rock 


a7 


CODE 


Z20 
Za 
Bae 
223 


1 Oy 


RIVERS. LAKES, AND CANALS 
Intermittent stream 
Intermittent river 
Disappearing stream 
Perennial stream 
Perennial river 
Small falls 
Small rapids 
Large falls 
Large rapids 
Masonry dam 
Dam with lock 
Dam carrying road 
Intermittent lake or pond 
Drv lake 
Narrow wash 
Wide wash 
Canal. flume. or aqueduct with lock 
Elevated aqueduct. flue. or conduit 
Aqueduct tunnel 
Water well 
Spring or seep 


GLACIERS AND PERMANENT SNOWFIELDS 
Contours and limits 
Form lines 


08 


CODE Deo ie ON 


Le) 
oO? 
= 


SUBMERGED AREAS AND BOGS 
Marsh or swamp 
Submerged marsh or swamp 
Wooded marsh or swamp 
Submerged wooded marsh or swamp 
Rice field 
Land subject to inundation 


Oo O 
N 


O 
Ho 


OW Wb lV bo 
oO o> 
Or co 


Oo 
Oo 


* Features not supported by the editor. 


Note : The editor does not provide automatic symbol plotting. Features containing curved lines 


can only be approximated using small line segments. 


begin UPDATE DATABASE algorithm 
{ 
initialize variables: 
initialize map region (TOTAL ROWS x TOTAL COLUMNS) 
pixel colors to unknown: 
initialize scratch pad values to 0: 
create a border around map region: 


/* index through each area of the map */ 
for (i = 0:1 << MAX DIMENSION; ++) 


{ 


new area = mapindex{i].next area: 


/* examine all areas at the same index reference point */ 
while (new area = exist) 


{ 


change new area colortype to new area description colorshade: 


/* increment to next area */ 
new area = new area -> next area: 


} 
j 


redraw map on screen: 
read color values within map off the screen: 
convert colorshades to appropriate processing color; 


call data compression algorithm: 


j 
end UPDATE DATABASE algorithm 
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begin DATA COMPRESSION algorithm 
{ 


initialize variables: 


/* create file to hold compressed data representation */ 
create output file: 


/* successively examine each pixel of the map region */ 
for (1 = 1:1 <= TOTAL ROWS; i++) 


{ 
for (j = 1:j <= TOTAL COLUMNS; j++) 


/* if pixel has not been examined as yet */ 

if (pixel(i, j) = not used) 
convert pixel color back to description colorshade: 
determine true map color: 
determine description of map area: 


write true map color to output file: 
write description code to output file: 
write speed code to output file: 


/* expand then subdivide homogeneous areas */ 
call expand area algorithm: 
call subset area algorithm: 


write end of area marker to output file: 


} 
j 
} 


j 
end DATA COMPRESSION algorithm 


Data Compression Algorithm 
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begin EXPAND AREA algorithm 
{ 


set pixel(x. y) = used: 
set corresponding (x, y) location in scratch pad to 1; 
insert pixel (x. y) location into the queue; 


/* examine only pixel values within the queue */ 
while (queue = not empty) 


{ 


remove (x, y) value off the queue; 


/* check the surrounding eight neighbors of the grid’s center to 
see if any are similiar in color tu 
for (i = (x-1): 1 <= (x41): i4+4) 
{ 
= for Gy yee 


{ 


/* tag map region and scratch pad location of similiarly 
colored neighbor. Make queue entry of pixel’s location */ 
if ((pixel color(i. j) = pixel color(x, y)) and (pixel(i, j) = not used)) 


set pixel(i. j) = used: 

set corresponding (x, y) location in scratch pad to 1; 

set Maximum x value = maximum(previous x value, 1); 
insert pixel (x, y) location into the queue; 


} 
} 
j 


end EXPAND AREA algorithm 


Expand Homogeneous Area Algorithm 
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begin SUBSET AREA algorithm 


{ 
do 


{ 


/* initialize temp variables to selected positions within the 
scratch pad vy 

set tempO = scratch pad area (x, y+1): 

set templ = scratch pad area (x+1, y+1): 

set temp2 = scratch pad area (x+1, y): 

set temp3 = scratch pad area (x+1, y -1): 


if (temp0 = 1) 
if (templ = 1) 


if (temp 2 = 1) 

{ 
expand type two polygon: fe VigtiresSegaitem b */ 
write triangle verticies to output file: 


j 


else 


{ 


expand type four polygon: /* Figure 5.3. itemd */ 
write triangle verticies to output file: 
} 
} 
else if (temp2 = 1) 
{ 


if (temp3 = 1) 


expand type five polygon; /* Figure 5.3, iteme */ 
write triangle verticies to output file; 


j 
} 


Algorithm for Subdividing an Area into Polygons 





else 


{ 
expand type six polygon: /* Figure 5.3, item f */ 
write triangle vertices to output file: 
} 
} 
else if (temp2 = 1) 
{ 
if (temp1l = 1) 


{ 


if (temp 3 = 1) 


expand type one polygon: /* Pigure oo ivenmieassey, 
write triangle verticies to output file; 


} 


else 


{ 


expand type three polygon: /* Figure 5.3. itemc */ 
write triangle verticies to output file: 
} 
} 
else if (temp3 = 1) 
{ 
expand type five polygon: /* Figure 5.3. iteme */ 
write triangle verticies to output file: 


} 


else 

{ | 
expand type eight polygon; /* Figure 5.3, item h */ 
write line endpoints to output file: 


j 
} 


Algorithm for Subdividing an Area into Polygons (contd) 
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else if (temp3 = 1) 

{ 
expand type nine polygon: } eee igureiorc. item 1) / 
write line endpoints to output file: 


} 

else if (templ = 1) 
expand type seven polygon; /* Figure 5.3, item g */ 
write line endpoints to output file: 


} 


else 


expand type ten polygon: / Se ihigure 5.3.atem | aaa 
write point location to output file: 


} 


~ reset scratch pad region flags: 


} 


while (x <= maximum x value); 


end SUBSET AREA algorithm 


Algorithm for Subdividing an Area into Polygons (cont'd) 
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